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C _~) Abstract. Laplace transforms for integrals of stochastic processes have been known in ana- 

£N| lytically closed form for just a handful of Markov processes: namely, the Ornstein-Uhlenbcck, 

I * the Cox-Ingerssol-Ross (CIR) process and the exponential of Brownian motion. In virtue of 

their analytical tractability, these processes are extensively used in modelling applications. In 
this paper, we construct broad extensions of these process classes. We show how the known 
models fit into a classification scheme for diffusion processes for which Laplace transforms 
OO for integrals of the diffusion processes and transitional probability densities can be evaluated 

as integrals of hypergeometric functions against the spectral measure for certain self-adjoint 
' , ' operators. We also extend this scheme to a class of finite-state Markov processes related to 

Ph hypergeometric polynomials in the discrete series of the Askey classification tree. 

Oh 

& 1. Introduction 

-t— > 
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where t < T, i? e C and (j), q : R — > R two Borcl functions. In this paper, we address the 
question of whether it is possible to compute the Laplace transform LT-t(Xt,i9) in analytically 
closed form. Our work builds upon several streams of research often motivated by applications 
to various fields of Physics and Finance, and unifies them to obtain a broad classification scheme 
for Laplace transforms expressible in analytic closed form. 
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Let (A" t )i> be a time-homogenous, real-valued Markov process on the filtered probability 
space (ft, {jF t } 4 >o,P) and consider the Laplace transform LT-t{X tl '&) defined as follows: 

(1.1) L T _ t {X u $) = E p f e -of t T 4>(x s )ds q{Xr) I Tt 



For q = 1 , a class of examples for which analytic closed form solutions are available is repre- 
sented by the so called affine models which are characterized by a representation of the form 

(1.2) L T _ t {X u §) = e ™(T-t,0)Xt+n(T-m _ 

The archetypical affine models are based on diffusion processes and are described by stochastic 
differential equations of the form 

(1.3) dX t = {a- bX t )dt + aX^dW t 

where a, b, a are constants and (3 — or \. The case (3 = corresponds to the Gaussian Ornstein- 
Uhlenbeck process (?) and the case (3 = \ corresponds to the Cox-Ingcrsoll-Ross (CIR) process 
(?). The case of the CIR process was generalized to bridges by Pitman and Yor in (?). It 
has been shown in (?), (?) and in (?) that any affine process which is a time-homogenous, 
nonnegative diffusion is necessarily of the CIR type. However, there are also affine processes 
with jumps. General non-negative affine processes correspond to the so called conservative CBI- 
processes (continuous state branching processes with immigration) and have been well studied, 
among others, by Kawazu and Watanabe in (?) and Filipovic in (?). 

An extension of the affine class, known as the quadratic class, postulates the Laplace transform 
being of the form 

(1.4) L T -t(X t J) = e l (T-t.V)Xt+m(T-t.V)X t +n(T-t^) ^ 
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The first examples of quadratic models appeared in the double square root model of LongstafF 
in (?) and in the nonlinear equilibrium model by Beaglcholc and Tenney in (?). Rogers (?) 
also uses examples where the pricing kernel is a quadratic function of the Markov process. Most 
recently, Filipovic (?) proved that if one represents the forward rate as a polynomial function of 
the diffusion process, the maximal consistent order of the polynomial is two. Consistency in this 
context means that the interest rate model will produce forward rate curves belonging to the 
parameterized family. Finally, Lcippold and Wu (?) formulated a general asset and derivative 
pricing framework for the quadratic class. 

A separate class of models for which the Laplace transform can be expressed in analytically 
closed form is represented by the exponential Brownian motion of equation 

(1.5) dX t = fiX t dt + oX t dW t 

where \i, a > are positive constants. This case was first considered by Yor in (?) who arrived to 
an expression involving a triple integral. An earlier related result for bond prices given in terms 
of an integral over modified Bessel functions was formulated by Dothan in (?). As an alterna- 
tive, Geman and Yor (?) derive a closed-form expression for the Laplace transform in terms 
of confluent hypcrgcomctric functions (see Donati-Martin et al. (?) and Yor (?) for further 
references). For applications to finance, one needs to compute the inverse Laplace transform, 
for which numerical methods have been developed by Geman and Eydeland (?), Fu ct al. (?), 
Craddock et al. (?), Shaw (?). Dufresne in (?) and Linetsky in (?) develop analytical methods 
and alternative expansions. 

In this article, we obtain far reaching extensions of the representation formula for the Laplace 
transform for the integrals of stochastic processes over geometric Brownian motions. The key 
idea is to seek expansions of similar form as those in (?) and (?) but expressed in terms of 
more general hypergeometric functions, see (?), as pioneered by Wong in (?). Let us recall that 
general hypergeometric functions are denoted as follows: 

(1.6) p F q (a 1 ,...,a p ;j 1 ,...,-f g ;z) 

for p < q+ l,7j G C \ — Z + , and are represented by the following Taylor expansion around z = 0: 

(ai)„...(a p )„ z n 



(1.7) pF g (ai,...,a p ;7i,...,7 g ;z) = ^ 



(7l)n ■ • • {lq)n Tl\ ' 

The Kummer functions in the work by Geman and Yor (?) are in the family \Fi of the so called 
confluent hypergeometric functions. Gaussian hypcrgcomctric functions arc in the family 2 F\ 
and admit the functions of type \F\ as limits. Both Gaussian and confluent hypergeometric 
functions solve differential equations of the Fuchsian class, see (?). Specifically 

(1.8) z(l - z) 2 F{'(a, 0; 7; z) + (7 - (1 + a + 0)z) 2 F{(a, /?; 7; z) - af3 2 F 1 (a, 0; 7; z) = 
and 

(1.9) z xFf (a; 7; z) + (7 - z) x F[(a- 7; z) - a ^(a; 7; z) = 0. 

In general, higher order hypergeometric functions are not associated to a differential equation. 
However, in some particularly important cases, they provide solutions of finite difference equa- 
tions. The Askey classification scheme, see (?) and (?), gives a complete list of all orthogonal 
polynomials solving either a differential or a finite difference equation and in addition satisfy a 
recurrence relation. All of these polynomials descend as particular or limiting cases from the 
so-called Racah polynomials, which are particular cases of the hypergeometric functions 4F3. 

We first consider the case of diffusion processes and next the case of finite state Markov pro- 
cesses. In the diffusion case, we construct a classification scheme based on reduction to eigenvalue 
problems admitting solutions within the class of Gaussian and confluent hypergeometric func- 
tions 2 F\ and \F\. In the second case, the problem is more difficult for several reasons, as there is 
no discrete equivalent of a theory of Fuchs type equations and, in addition, the groups of confor- 
mal transformations and diffcomorphisms do not extend to lattices. What we do in the discrete 
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case is to take the moves from the Askey classification scheme for orthogonal polynomials and 
show how to extend the previous spectral decomposition to the case of Meixner, dual Hahn and 
Racah polynomials, which are special cases of 2-F1, 3F2 and 4-F3 hypergeometric functions. 

For a diffusion process, on a domain D x C M, of the form 
(1.10) dX t = fi(X t )dt + a{X t )dW u 

we define the transitional probability density pT-t{x, y) as the density of the Markov semigroup 
of the process X t : 

(l.H) E p [ f(X T ) \T t ] = [ f(y) p T _ t (x,y)dy. 

J D 

We arc interested in building a classification scheme for the drift and volatility functions fi(x),a(x) 
such that the calculation of both functions p and L can be reduced to computing an integral 
over hypergeometric functions. The transitional probability density and Laplace transform can 
be computed in terms of the spectral resolution for the infinitesimal generator of the process X t 

^^Tdx-^^dx 

and the Feynman-Kac operator 

(1-13) £ = = ^~ + K*)~- H{x). 

In fact, we have that 

(1.14) p T _ t (x,y)=e^ j: (x, y ) 
and 

(1.15) L T _ t {x^)= [ q(y) e^ T - t ^(x,y)dy. 

J D 

As we show in detail in Section[2j these operators are conjugated by a non-singular transformation 
to self-adjoint operators which admit a spectral resolution. The calculation of the transitional 
probability density in ( 1.14| and the Laplace transform in ( 1.15 1 is thus reduced to the resolution 



of the differential eigenvalue problems 
(1.16) Cf(x) = Xf(x) , £f(x) = \f(x). 

To properly define the classification scheme, we specify by what means the reduction can be 
accomplished. 

Definition 1. The problem of finding the transitional probability density and the Laplace trans- 



form for the process in (1.10) is said to be reducible to a spectral integral over hypergeometric 



functions if the two eigenvalue problems in {1.16) can be recast in the form of a differential 



equation for hypergeometric functions such as either (1.9) or {1.8), by means of a combination 
of the following three operations TV 

(1) Tz- change of variable x t— > z = Z(x) where Z(x) is a diffeomorphism Z : D x — > D z 
such that 

(1.17) C X ^C Z , C X ^C Z , 

(2) Th'. gauge transformation associated to a strictly positive function h such that 

(1.18) C^h^Ch , C^h- l Ch, 

(3) T 7 2 : left-multiplication by a strictly positive function 7 2 such that 

(1.19) £^ 7 2 £ , £(-> 7 2 £. 

The third kind of transformations was first recognized in its generality by Natanzon in the 
article (?) on integrable Schrodinger equations, see also Milson's paper (?). The following 
theorem gives a concise statement of our main classification result: 
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Theorem 2 (First Classification Theorem). The most general reducible diffusion process (up to 
diffeomorphism) according to Definition^ can be constructed as follows: 

(1) four second order polynomials in x: A(x), Q(x, R(x), S(x), such that A(x) belongs 
to the set {1, x, x(l — x), x 2 + 1} and R(x) > 0; 

(2) conditions for the stochastic process X t on the boundary of the domain D x , specifying 
the relative probability of reflection versus absorption upon hitting the boundary; 

(3) a solution of the following equation in D x for some £ G K: 

(1.20) A j ^h"(x) + §?\h(x)=th(x) 

R(x) R(x) 

The function h(x) is a linear combination of hypergeometric functions of the confluent 
type i Fx if A{x) G {l,a:} and of the Gaussian type %Fx if A{x) G {x(l — x),x 2 + 1}. 

The process associated to this choice is given by a generic solution to the following stochastic 
differential equation on the domain D x : 

with the boundary conditions above. The Laplace transform is specified by 

(L22) ^ {x) =mx)- 



The proof of this theorem is in Section [2] These constructs are based on spectral analysis 
techniques for which we refer to the book by Reed and Simon (?). We also make the spectral 
analysis more explicit and list the expressions for the transitional probability density and Laplace 
transforms in terms of the kernel of semigroups generated by integrable quantum Schrodinger 
operators. In Section [3] we specialize further and re-discover the known cases of processes built 
upon the geometric Brownian motion and on the Ornstein-Uhlenbeck and CIR processes, along 
with some interesting extensions. 

In Sections [4] and [5] we restrict the framework to the special case of hypergeometric polyno- 
mials. In the discrete case, studied in Section [5] the process X t takes on only a discrete set of 
values, as opposed to following a diffusion process. In this class of models, we base our analysis 
on the Askey- Wilson theory of orthogonal polynomials, see (?) and (?). We briefly review the 
basic notions, following (?). 

Definition 3. An orthogonal system of polynomials is given by a sequence of polynomials Q n (x) 
of order n for n G N on the interval D C K which satisfies an orthogonality condition of the 
form 

(1.23) / Q n (x)Q m (x)p(dx) = d 2 n 5 nm , n,meN, 
Jd 

where the d n are constants and p(dx) is a given measure. One distinguishes between continuous 
polynomials whereby p{dx) is absolutely continuous with respect to the Lebesgue measure, i.e. 

(1.24) p(dx) = w(x)dx 

for some weight function w{x), and discrete polynomials for which 

N 

(1.25) p(dx)=^2w(xi)S{x-i) , N eN. 

i=0 

All orthogonal polynomials satisfy a three-term recurrence relation of the form 

(1.26) xQ n (x) = A n Q n+1 (x) - B n Q n (x) + C n Q n ^x(x) 

where n > 1, A n > 0, C n > and B n G K. Together with the conditions Q_i(x) = and 
Qo(x) = 1, all the Q n {x) can be determined based on this recurrence relation. The converse 
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is also true and is known as the Favard theorem, see (?). Moreover, they satisfy the following 
eigenvalue equation, 

(1.27) CQ n {x) = X n Q n (x), 

for C a second-order differential operator in the continuous case or a finite difference operator 
in the discrete case. 

The reducibility condition in Definition[T|is mirrored by the following (inequivalent) one which 
refers to orthogonal polynomials in the continuous series as opposed to Gaussian hypergeometric 
functions: 



Definition 4. The problem of finding the transitional probability density in (1.14) and the 



Laplace transform in (1.15) is said to be reducible to a spectral integral over orthogonal poly- 



nomials if the two eigenvalue problems in {1.16) have the same orthogonal polynomials as eigen- 
functions. 

We first restrict the framework to continuous orthogonal polynomials that have as generator 
a second-order differential operator, 

a 2 d 2 d 

(1.28) C -Y A ^d X 2+ ^~ bx ^ 

acting on the Hilbert space L 2 (D,p), where A(x) S {l,x,x(l — x)}, a € R and b, a > 0. This 
class consists of the Hermite, Laguerre and Jacobi polynomials up to diffeomorphism. Our main 
result concerning the latter continuous orthogonal polynomials can be stated as follows: 

Theorem 5 (Second Classification Theorem). The most general reducible diffusion process (up 
to diffeomorphism) in the sense of Definition^ has infinitesimal generator C given by (1.28). 
Its transitional probability density can be expressed as 

°^ e A„(T-t) 

(1.29) Pr-t(x, y) = — -p Q n (x; a, b)Q n (y; a, b)w(y). 

n=0 " 

Furthermore, for some parameters a, b, C S R, 

(1. 3 0) - c + «. - to) - , 5 - b,)^L + (5 -^- fa) ' 

and the Laplace transform is given by the following convergent series: 

(a-by)-(a-by) ^ j^^T-t) 

The coefficients z n are given by: 

(1.32) Zn = ^J D q ^ exp (- f {a ~ b a2A(y)~ bV) dy ) Qn{x '' ^ l)P{dx) - 

In Section [4] we present the proof of this alternative classification scheme based on orthogonal 
polynomials of the continuous series. This discussion sets the premise for the extension of the 
result to orthogonal polynomials in the discrete series. Discrete orthogonal polynomials are 
characterized by a finite difference generator on the Hilbert space 1 2 (An,w) with Ajv the set 
{0,...,N}. 

Definition 6. Let A h and denote the difference operators defined as follows: 

(1.33) A h y(x) = y(x + h) - 2y(x) + y(x - h) , V\y{x) = y(x + h) - y{x). 
Using these operators, the finite difference generator C takes the form 

(1.34) C = -D(a;)A 1 + (D(x) - B(x))V+, 



(1.31) L r _ t (s,#)=exp( / 1 J 2 A( y) dy l J2 eXniT ~ t)z M^a,b)- 
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where B(x) and D(x) are rational functions of at most fourth order in the numerator and at 
most second order in the denominator. Our main result in the discrete case can be stated as 
follows: 

Theorem 7 (Third Classification Theorem). The most general reducible discrete Markov process 
in the sense of Definition^ has infinitesimal generator C given by (1.3^). Its transitional 
probability density can be expressed as 

A e A -( T -*) 

(1-35) pr-t(x,y) = — — 2 Qn(x)Q n (y)w(y). 

71=0 " 

4>{x) must be of the form 

(1.36) 4>{x) = B(x) + D(x) - B(x) - D(x), 

where B(x) and D(x) are the same rational functions as B(x) and D(x) up to a multiplicative 
constant, but for different parameters, and satisfy the lattice condition 

(1.37) B(x- l)D(x) = B(x- l)D(x). 

o 

The Laplace transform is given by the convergent series — 1 by convention): 

k=l 

(1-38) L T _ t (z, 1) = f[ §| E ^^ZnQnix). 

fc=l ' n=0 

The coefficients z n are as follows: 

(1-39) z n = i f[ rM l{*)Qn{x)w{x). 

n x£A N k=l ^ ' 



In our notation, {Q n (x)} is the same set of orthogonal polynomials as {Q n (x)} expect for the 
value of the parameters. Section [5] gives a proof of the latter result, as well as explicit represen- 
tations for processes based on the Meixner, the dual Hahn and the Racah polynomials. 

This paper is organized as follows: Section [2] presents the proof of the first classification 
scheme in the diffusion case and goes in more detail to provide a spectral representation formula 
for transitional probability densities and the Laplace transform. Section [3] contains a discussion 
of the classical examples showing how the geometric Brownian motion, the Ornstcin-Uhlenbeck 
and the CIR process fit in this classification scheme. Section [4] contains the proof of an alterna- 
tive classification scheme based on orthogonal polynomials of the continuous series, along with 
interesting examples. This discussion sets the premise for Section [5] where we extend the result 
to orthogonal polynomials in the discrete series. Finally, in Section [6] we discuss limiting rela- 
tion and establish connections between models corresponding to the discrete and the continuous 
series, which is a useful result to construct numerical very stable discretization schemes. 
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2. Classification Theorem for Diffusion Processes 

In this section, we prove our first classification result, Theorem [2] in the diffusion case. We 
start by reviewing some background notions concerning Fuchsian differential equations and the 
so-called Bose invariants, and then proceed to the proof of the theorem. 

2.1. Fuchsian Differential Equations. Consider the second order partial differential equa- 
tions for the holomorphic function F(z) 

(2.1) F"{z) + p{z)F'{z) + q(z)F(z) = 
for some holomorphic functions p(z),q(z). 

Definition 8. Let a £ C be an isolated singularity for the holomorphic function F{z). The 
singularity in a is called regular if there is an exponent p £ C for which the function (z~a)~ p F(z) 
admits a Laurent expansion with finitely many negative powers around z — a, i.e. 

o 

(2.2) F(z) = (z-ay J2 (*"«)" 



n— — m 



for some m G N. The point oo is a regular singularity of the function F{z) if z = is a regular 
singularity of the function F (-) . 

In (?), Fuchs gives conditions on the coefficients p(z) and q(z) which ensure that solutions 
have only regular singularities. 



Theorem 9 (Fuchs). Let F(z) be a solution of equation (2.1) with singularities in the points 
a\, . . . ,a n and oo. Then these singularities are all regular if and only if the functions p{z) and 
q{z) have the form 

( 2 -3) P{Z) = r r 

{z-a 1 )...(z- On) 

and 

( 2 - 4 ) = 7 ^ « 

(z - ai) z ...{z- a n y 

where Po(z) is a polynomial of order (n — 1) and qo(z) is a polynomial of order 2n — 2. 

An alternative expression for the coefficient p{z) of an equation with only regular singularities 

is 

Si 



(2-5) P (z) = 



Z-Cti 



6XP Q j p ( w ^ dw ^ = CjJ^(z - OiY 



where the Si, i = 1, ..n, are constants. In particular, we have that 
(2.6) 

where C is a constant. The function 

n 

(2.7) F(y) = l[(y-a l )^F(y) 

i=i 

solves the equation 

(2.8) F"(y) + I(y)F(y) = 



where 



(2-9) I(y) = --p'(y) + -p( y y + q(y). 
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Definition 10. The function I{y) is called the Bose invariant of the equation (2.1). Notice that 
I(y) has the form 

(2-10) i(y) = , ^ V2 

(y-a 1 y...{y-a n y 
where Io(y) is a polynomial of order In — 2, without restrictions on the coefficients. 

Let us focus again on the case n — 2, assume that coefficients are real and that only real 
linear fractional transformations are allowed to move the regular singularities. In this situation 
we have to distinguish between two different cases for Bose invariants: 

• Case I (ai = 0, a 2 = 1): 

(2 ' n) 1{v) = Ti^W 2 

• Case II (ai = i, a 2 = —i): 

, 010 s, r ,v so + siy + s 2 y 2 

(2 ' 12) 1{y) = W + 1) 2 



These cases reduce to the Gaussian hypergeometric equation (2.18 1 for the function 2 Fi as is 
shown below. Furthermore, special cases for the Bose invariant occur in the limit when either 
ct\ or a 2 or both roots tend to oo, i.e. 

• Case III (a\ = 0, a 2 — oo): 

(2.13) /(,) = S0 + ^ 2 W 

y 

• Case IV (ai = oo, a 2 — oo): 

(2.14) I(y) = s + s iy + s 2 y 2 

Case III reduces to the confluent hypergeometric equation iFi and Case IV corresponds to the 
case of triple confluence at infinity. Notice that the above four cases can all be captured by a 
single expression as stated in the following: 

Remark 11. The Bose invariants corresponding to the Gaussian hypergeometric function 2 F\ 
and to its confluent limit can be reduced to the following normal form by means of a real valued 
linear fractional transformation: 

(2.15) m - H 

where A(y) £ {y(l — y), y 2 + 1, y, 1} and Q is a polynomial in y with deg Q < 2. 

The first two cases correspond to three regular singularities at distinct points. In these cases, 
solutions can be expressed through Gaussian hypergeometric functions 2 Fi. Fractional linear 
transformations of the form 

az + b 

(2.16) z ~ — - - 

cz + a 

where a,b,c,d£ C and ad — be ^ 0, are one-to-one maps of the extended complex line C U oo 
into itself and map regular singularities into regular singularities. By applying a fractional linear 
transformation, one can map the singularities ct\ and a 2 to and 1, respectively. Furthermore, 
we have transformations of the form 

(2.17) F(z) ^ (z - ai) Pl (z - a 2 ) P2 F(z). 

The combination of these two transformations allows one to reduce any Fuchsian differential 
equation with three regular singular points to the form 

(2.18) z(l - z)F"{z) + (7 - (1 + a + P)z)F'(z) - a/3F(z) = 0. 

The function 2 Fi(a, z) is an elementary solution of this equation along with 2 Fi(a, (3; 1 + 
a + [3 - 7;1 — z). 
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Case III corresponds to the limit when a regular singularity merges with the regular singularity 
at oo while the other one stays at 0. This limit can be obtained starting from the equation 
corresponding to two coinciding singularities at 0, i.e. ot\ = a 2 — 0: 

;"2 



(2.19) 



F"{z) + 



ci 



c 2 z_ F , {z) + c 3 + c 4 z + c 5 z F[z) = o 



By applying the coordinate transformation z 1— > - we find that 

(2.20) z 2 F'\z) + ((2 + c 2 )z + c lZ 2 )F'(z) + (c 3 + c 4 z + c 5 z 2 )F(z) = 

By rescaling the independent variable z and rescaling the function so that F{z) 1— > e pz f(ujz) this 
equation reduces to the Kummcr differential equation 

(2.21) zF"(z) + (7 - z)F'(z) - aF{z) = 

which admits 1^1(0;; 7; z ) as a solution. In alternative, one can reduce equation (2.201 to the 
form 

,2 < 



(2.22) 



F"(z) + 



A 

z 



F(z) = 



which is called the Whittaker differential equation. The case where all three singularities merge 
at 00 is also interesting and can be solved by rescaled confluent hypergeometric functions. 

2.2. Proof of the First Classification Theorem [2} We start by presenting obvious facts 
about the transformations Ti in Definition [T] 



Th 1 = T ± 



T~2 = T 1 
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Remark 12. The transformations Ti are invertible with respective inverse: 
(2.23) T z 1 = T x 

where X : D z — > D x is the inverse of Z{x). 

Remark 13. The transformations Ti commute with one another. 

The proof of the First Classification Theorem [2] for diffusion processes follows. The process 
X t has infinitesimal generator 



(2.24) 



C = 



A(x) 2 d 2 h'(x) A(x) 2 d 



R(x) dx 2 h{x) R{x) dx 

To solve the eigenvalue problem Cf = A/, define the following left-multiplication T 7 2 and gauge 
transformation T^: 



(2.25) 



T^C 



and T h C = \ Ch ' 

R(x) h 



so that the eigenvalue equation then transforms into 

d 2 S(x) 



(2.26) 



T i T -X 



dx 2 



A(x) 2 , 

As R(x), S(x) are second order polynomials and A(x) € {1, x, x(l — x), 1 + x 2 }, the solution f(x) 
can be expressed as a hypergeometric function with Bose invariant 

(2.27) I(x) = A{x)2 . 

The same operations can be applied to the eigenvalue problem for the Feynman-Kac operator, 
Cf = Xf with 

(2.28) C = C-&<t>(x), 
which leads to 

8 2 , S(x) 
dx 2 



(2.29) 



Cf=~Xf 



A{xf 



»,> { .n l/(x) = A^/». 



A{xf 
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For f(x) to be expressed as a hypergeometric function, we require 
(2-30) = jj^ 

for Q(x, i?) a second order polynomial in x. The Bose invariant for the equation is 

(2.31) I(x) = . 

The converse follows from the facts that the transformations are invcrtiblc and the Bose invari- 
ants are both expressed in the most general form. 

Remark 14. The first transformation (Ty: change of variable) has not been used in the proof. 
X t is therefore the most general reducible diffusion process only up to diffeomorphism. 

2.3. Spectral Resolutions. Theorem [2] describes all processes with explicitly solvable transi- 
tional probability density and Laplace transform for the integral of the process. We wish now to 
give a closed form expression for both quantities. The following lemma allows one to determine 
the nature of the spectrum of the operators C and C. The nature of the spectrum is indeed based 
on the shape of the Schrodinger potential, coming out of the eigenvalue problem once reduced 
to a Schrodinger equation. Theorem [16] shows how this transformation operates on the kernel 
of the semigroup generated by the Schrodinger operators and gives a general spectral resolution 
of the operators C and C. 

Lemma 15. Let T — T„ Tz T7 1 where the diffeomorphism Z : D r — > D z is qiven by Z'(x) = 

A(x) V R( x ( z )> J 

operators C and C reduce to the following Schrodinger operators: 

TC = ~ - U x {z) = -Hi 
oz z 

(2.32) TC = ^-U 2 (z) = -U 2 



dz 2 

g'V fg'\ S(X(z)) 



gj \gj R(X(z)) 
g'\ 2 /g'V S{X(z))-Q(X(z),d) 



where the potentials are given by 
U(z) = 

(2-33) U(z) = , 

\9 ) \9 ) R{X{z)) 

and ' denotes the derivative with respect to z. 
Proof. 

T g T z T h C - Tg T z y 1 ^^ + 1 ^- ) 

T Z " 9 | M 

9 \dz 2 {Z'fdz R(X(z)) 

dz 2 g 

and 



dz 2 g \gj R(X(z)) 
T g T z T- X L = T g T z T-'C - ^(X(z)). 



a 



The two Schrodinger operators Hi and H2 defined in the previous lemma have a spectral 
resolution, 

(2.34) H^ p (z) = p $ p (z), 
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given by a complete set of normalized cigcnfunctions <f> p (z) for p = —A or p = —A, z = 1,2 
respectively. The spectrum is in general a combination of a pure-point spectrum cr pp (M^ ) and an 
absolutely continuous spectrum cr ac (IHIi). The kernel of the semigroup generated by the respective 
Schrodinger operators has the general form 
(2.35) 

e -(T- t )H l(zf);Zi)= e-( T - t ^ p (z )^ p (z 1 )+ f e-V-^riz^^dkip) 

with dk(p) — — , dp and U~ the lowest limit of the potential Ui(z) as z tends to the boundaries 

2V P~U i 

of the domain D z . 

Theorem 16. The transitional probability density and the Laplace transform of any reducible 
process described in Theorem^ by the operator C and L is related to the kernels of the semigroups 
generated by the respective Schrodinger operators as follows: 

^2\ 1/4 / Df „ A \ 3/4 



^' K <-> - m m) (H) ^°" 1 <*«' 

,, 36 , n,y, - ^ (^) V4 (™) V4 ^(* W ,z M ). 

Proof. The transformation T defined in Lemma [15] extends to the kernels as follows: 

and similarly for C and H2. Recall that g(Z(x)) = I — — — ) and the Jacobian of Tz is given 

V R \ x ) J 

by — = , which concludes the proof. □ 

dy A(y) 
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3. Examples of Solvable Diffusions 

In this section, we show that the transitional probability density and the Laplace transform 
for the integral of the geometric Brownian motion, the Ornstein-Uhlenbeck process and the CIR 
process arise as corollaries of the First Classification Theorem [2] 

3.1. The geometric Brownian motion. 

Definition 17. The geometric Brownian motion is defined by the solution of the following 
stochastic differential equation: 

(3.1) dX t = fJ,X t dt + aX t dW t 

with initial condition X t= o = xq. 

Corollary 18. The transitional probability density for the geometric Brownian motion is given 
by the following formula: 

1 x 



(ln(|)-( M -4)(T-f))' 



(3.2) p T _ t (x,y) = — \ exp 

y^2ira 2 {T-t) I 2a\T-t) 

The Laplace transform is explicitly solvable if and only if 



(3.3) ^( !B ) = _^^i_^j_ to _ tlX + t2a; 
where to,t\ € K and t% > could depend on #. It is then given in terms of the Laguerre 

N 

polynomials L„ and the Whittaker function (by convention, = if N < 0): 

n=0 

N 

Lr-t(x,0) = x-$ + i e~^* J2 z " e(T ~ t)Kx ^ L n n) ( 2 ^ x ) 

n=0 

(3.4) + x-5* \ z fc e- (T - t)(fc2+t; - ) M h s,(2Vhx)dk 

JO 2^*2' 2 

where N = — g] denotes the integer part oft), S n = ~2n — 1 + U— = — to) 
and (5fc = iy^zk. For n = 0, 1, . . . , JV, £/ie discrete eigenvalues are given by 

^ A ^y((" + 2-2^) 2+t °-4)- 



T/ie coefficients z n and Zk are respectively 

T(<5n) io 



(3.6) z„ = f%(x) ^ + "^^4^(2^)^ 



and 



1/2 /" 

(3.7) z k=^\^T q{x)x^~ 2 M H s k (2Vhx)dx. 

2vr V cHt 2 Jo ^7*?' 2 

Proof. The infinitesimal generator is 

„ ff 2 2 9 2 3 
(3-8) C =Y X dxi+^dx 

whereas the Feynman-Kac operator has the form 

(3.9) £=^^ +M *A +#(,). 
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This case fits the classification scheme in Theorem |2] if one selects 

2 h! p 



(3.10) 



A(x) = x , R(x) = 



This choice sets the shape of the polynomial S(x), from (1.201 in Theorem|2] to be 



(3.11) 



S(x) 



1 



2Z + ft(l 



which in turn defines the Bose invariant as 

(3.12) /<*) . + -,('-*> . 
The Schrodinger potential, given by 

(3.13) U(z) = - 



,('-£)• 

is constant for all z. Hence the spectrum is absolutely continuous. The normalized eigenfunctions 
for the Schrodinger equations (2.34) are 

(3.14) $,(*) 



1 



/2n 



where k 2 = p — ^- 
the operator C: 



(3.15) 



f (1 



(x,y) 



Theorem 



1G 



yields the kernel for the semigroup generated by 
( m (2)-( M _4)(T-t)) ; 



^/2Tia 2 {T-t) 



exp 



y 



2a 2 {T-t) 



which is pT-t{x,y) and in which one recognizes the transitional probability density of the geo- 
metric Brownian motion. 

The kernel of the semigroup generated by C is more general, since there is no restrictions on 
the polynomial Q(x,"&), which, for to,t%,t2 € K, can be written in the form 

(3.16) Q(x,d) = S(x) - t - hx+ t 2 x 2 . 

The Bose invariant in this case is 

2A 



(3-17) 



I(x) 



-t 2 + ~ 



to 



which gives rise to two independent solutions to the eigenvalue problem u" + Iu = 0: 
u±(x) 



(3.18) 



where -j- 



(3.19) 



M h ± so (2Vhx) 

2^*2' 2 
, ±5n + l it- 

(2Vhx)^e~^ x 1 F 1 



±s +i _ ti 

2 

±<y + 1 



2y / t 2 ~X 



to + ^r. The Schrodinger potential, given by 



U(z) = 



1 



- to - Ue^ 



+ t 2 e 



has no singularities and is bounded from below if t 2 > 0. It indicates that the spectrum is not 
strictly discrete, since U(z) ^ooasz->oo but U(z) — »• Z7_ e 
the description of the spectrum separates in two cases: 

• If ti < 0, then U(z) is monotonously increasing on D z 
uous: p = —A > [/_. 



— to) as z — » — oo. Hence, 
= K and the spectrum is contin- 



• If ti > 0, then U{z) has a minimum, 



%rirhi at z — \ / In tjt" and the 

2 412 V °" 2 



spectrum is discrete for Uq < p < U- and continuous for p > U- 
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Both cases can however be solved simultaneously. The solution to the Schrodinger equation, 
$p(z), is related to u±(x), via the diffeomorphism x = X(z) = eV^ z and the gauge transfor- 
mation g{z) = (^-)ieV / ^ z , as follows: 

(3.20) <S> p {z)=g- 1 {z)u ± {X{z)). 

The asymptotic behavior as z — > — oo, $ P {z) ~ ~^ e±lkz with k = \fp — U- > 0, enforces the 
normalization condition for <& p {z) to the following: 

(3.21) = ^ r (2Vt 2 )- 1 -^e-^ z u ± (X(z)). 

The continuous spectrum appears for p greater than the lowest limit J7_, implying that Sk = 
5o(k) — i\J~^k is imaginary. Whether or not the discrete part of the spectrum has an infinite 
number of discrete levels depends on the asymptotic behavior of the potential. As z — > — oo, 
x — ► and U(Z(x)) develops up to second order: 

2 / i 

(3.22) U(Z(x)) = — (--t -t 1 x + t 2 x 2 

Since t\ > for the discrete spectrum, the convergence to the limit XJ_ as x — > is faster than 
x 2 , which implies that there is only a finite number of bound states. 

The hypergeometric functions in the solutions u±(x) reduce to polynomials for x <E [0,oo) if 
respectively 

(3.23) + 1 - -4= = -„ e -N. 
V ; 2 2JT 2 



We set S n = So(n) to emphasize the dependance on n. The solution u-(x) is not L 2 -normalizable, 
whereas m+ (x) is given in terms of the Laguerre polynomials 



(3.24) u n (x) = C{2^T 2 x) S - a ^ 1 e~" /r2X L^\2^r 2 x) 

/oo 
|<I> rl (z)| 2 <iz = 1 fixes the 
-oo 



,-2 



constant C = * — — - (cf. (?), p. 462). The discrete eigenvalues p = — A„ are given by 

y 4Vt 2 r(d„) 

(3.5 1, where n is restricted to {0, 1, ... , N} by the condition — A„ > Uq, and N = |"-^= — ~] . 

Hence, from Theorem 16 the kernel of the semigroup generated by C is given by the following 
spectral resolution: 

(3.25) eF-Qtfay) = aT^+iy^-i e -V*2(*+y) 

■Ee (r - t)A »^(^)-4 i " ) (2^)e>(2^) 

*— ' r(<5„) 



n=0 



1 / 2 



2tt V oH 2 



e -(T-t)(* fi , fc {2^/T 2 x)M t A - hL (2Vt 2 ~y)dk. 

Finally the Laplace transform of the latter kernel is given by integration over _D X = [0, oo) and 
yields (JO). □ 
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3.2. The Ornstein-Uhlenbeck process. In this subsection, we restrict the framework to the 
affine models, i.e. we set 

(3.26) q(x) — exp (u)<f)(x)) 

for some We show that in the special case of the Ornstein-Uhlenbeck process, both the 

transitional probability density and the Laplace transform can be expressed as summations over 
Hermite polynomials. 

Definition 19. The Ornstein-Uhlenbeck process is defined by the solution of the following sto- 
chastic differential equation: 

(3.27) dX t = (a- bX t )dt + adW t 
with b > and initial condition X t= o = xq. 

Corollary 20. The transitional probability density is given by the following formula: 

, 2 



(3.28) 



where z{x) = 



PT _ t(x , y)= J_(l 



,-2b(T-t) 



-1/2 



exp 



(z(y) - z{x)e-^ T -^)' 

1 _ e ~2b(T-t) 



The Laplace transform is explicitly solvable if and only if 



(3.29) 



4>(x) 



2d \(T 



ab — ab 



with to, a, b e K and could depend on It is of the quadratic form 

(3.30) L T _ t {x, 0, W) = e m{T-t)-n(T-t) X -l(T-t) X i 

where the functions rn(r), n(r) and l(r) are as follows 



1 



m ( T ) = o ln 



- I In ( P - (p - 



-2br 



a 2 b 



b- 



- <T Z t ( 



+ 



20 



o- h )+p[q + 



2pq^e- b - + (0- +pq 2 ^) e 



-2br 



p-{p~^)e 



b_\ c -2bT 



n(r) 

l(r) 
with 



a- a 2pq£e-» T + (p-£)2^e 



a n -2br 



b-b 
~2o^ 

P : 



b_\ -2bT 



(P" 



b_ p -2br 



P- (p 



(b+b)(oj(b-b) + $) 

2(T 2 $ 



<1 : 



??(a + a) + uj(ab — ab) 



{b + b)(&-u(b-b)) 
Proof. The infinitesimal generator of the Ornstein-Uhlenbeck process is 



(3.31) 



C= Ydx^ + {a ~ bx) dx- 



and the Feynman-Kac operator has the form 

- a 2 d 2 



(3.32) 



d 



2Qx 2+{a - hx) Qx-^- 



This case fits the classification scheme in Theorem |2] if one selects 

(3.33) A(x) = l , R(x)=% , h(x) = e^ x -^ x \ 

o~ 

This choice implies that the shape of the polynomial S(x), from (1.201 in Theorem [2j must be 

1 , a\ 2ab b 2 



(3.34) 



S(x) = ^(2£ + b- 



r) + 



-x 2 
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which in turn defines the Bosc invariant as 

1 



(3.35) I(x) 
The Schrodinger potential, given by 



(2£-2A + 6--) 



2ab 

V 4 



(3.36) U ^=-t-\ + &> 



ab b 2 , 

/2a 2 4 



with z = yf-^x, goes to oo as z — > ±oo. Hence, the spectrum displays discrete eigenvalues of 
the form A„ = —bn for n G N. The corresponding normalized eigenfunctions that satisfy the 



Schrodinger equation (2.34| are 
(3.37) $„(Z(x)) 



1 



Vn\2 n \2tt 



Theorem [16] yields the kernel for the semigroup generated by £ as a summation over Hermite 
polynomials. Using Mehler's formula (cf. (?), p. 710) and the notation z(x) = 



the kernel re-sums into 



(3.38) 



,{T-t)C 



i - e - 2b ( T -*: 



-1/2 



exp 



(z(y) - z{x)e- b ( T -V) 



21 



1 



,-26(T-t) 



which is pT-t{x, y) and in which one recognizes the probability density of the Ornstein-Uhlenbeck 
process. 

For convenience and without loss of generality, we set the form of the polynomial Q(x, i?) to 

„, „, „, , 2ab b 2 ~ 

(3.39) 



Q(x, ■d) = S(x) —t ^x + —x*. 



2ab 

cr 4 a 

where to,a,b could depend on The Bose invariant in this case is 

2A 2ab b 2 



(3.40) 



I(x) =t Q -— + 



The Schrodinger potential, given by 



(3.41) 



ab a 



U(z) = -— - —\i —z + 



is very similar to U(z) and indicates that the spectrum is again discrete since U(z) -t oo as 
z — > ±oo. The solution & n (y) is given in terms of the Hermite polynomials 



(3.42) 



Vn\2 n \ 27r 



1/4 



r(-f) : 



H„ 



for the eigenvalues A„ = — bn — | + ^ 2 , 
generated by C is given by the following spectral resolution 



%rto. From Theorem 



16 



the kernel of the semigroup 



,(T-t)C 



(x,y) 



exp 



a 2 b a- 



(a + a) 



(3.43) 



^° P A»(r-t) 



n=0 



2 n n\ 



b , <i . 
(x 



JL_ ( J + 6 ) + - (a _„, 

i^f^fr-!)V 



2d 



The integration of the latter yields the Laplace transform for q(x) = exp(wx) as a convergent 
series in terms of the Hermite polynomials which re-sums to the formula (3.301 (cf. (?), p. 
488(16) and p. 710(1)). □ 
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Remark 21. Computing the Laplace transform L in the case where 4>{x) is affine, i.e. 4>(x) = x, 
is a direct consequence of the previous corollary. Setting the parameters to 

b / a \ 2 a 2 -& — 

(3.44) f = - , a = a-— , b = b 



proves the following proposition: 

Remark 22. The Laplace transform for the affine Ornstein-Uhlenbeck process is as follows: 

(3.45) L T _ t (x, 1?, CJ) = e m{T-t)-n{T-t)x 

where 

d- {§ + ujb)e- bT 



(3.46) 



n(r) 

(n(r)+u;-i?r)( a b-^) a 2 2 2 
m(r) = p 4ftK T ) ~ w )• 



3.3. The CIR process. In this subsection, we focuss again on the affine models, i.e. we set 

(3.47) q(x) = exp(iox) 

for some wel. We show how to derive from the First Classification Theorem [2j the transitional 
probability density for the CIR process and the Laplace transform in the affine case. The main 
tool is the use of the Laguerre polynomials as cigcnfunctions for both the infinitesimal generator 
and the Feynman-Kac operator. 

Definition 23. The CTR process is defined as the solution of the following stochastic differential 
equation on D x = K + : 

(3.48) dX t = (a- bX t )dt + a^X~ t dW t 
with a, b > and initial condition X t= o — Xq. 

Corollary 24. The transitional probability density of the CLR process is given in terms of the 
modified Bessel function I a as follows: 



ye 



exp 



c I y + xe 



" b(T - 4) )] fey/xye-W-*)) 



Pr-t{x,y) = 
(3.49) 

with c = c(T — t) = — e~ h ( T ~ t ^)~ 1 . The Laplace transform is computable in closed form if 
and only if 

(3.50) <p(x) ■■ 



1 



a 2 t \ 1 ( ab a 2 ti 



<J 2 t 2 



2a 2 



x, 



where to, ti, ^2 € K could depend on i9. 

In the affine case, where 4>{x) — x, the Laplace transform is of the closed form 



(3.51) 

where 



(3.52) 



to(t) 
n(r) 



2a 



In 



L T ^ t (x,ti,u J ) = e m ( T - t ^ T - t ^ 

l e br/2 

ocosh(^) + (6 - uju 2 ) sinh(^) 



(b - loo- 2 ) 



2\2 



sinh(^) 



5cosh(f ) + (b - loo 2 ) sinh(f ) 



and b = V^a 2 + b 2 . 
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Proof. The infinitesimal generator of the CIR process is 
(3.53) C 



j 2 d 2 . u ^ d 
2 X d^ + {a - bx) di 



and the Feynman-Kac operator has the form 
(3.54) t : 



a- 8 2 . , , d 
i— + (a — ox) - 



■$<t>(x). 



2 dx 2 ' v ' dx 
This case fits the classification scheme in Theorem [2] if one selects 

2x 

O 2 

This choice sets the form of the polynomial S(x), from (1.201 in Theorem [2] to be 



(3.55) 



A(x) = 



R(x) 



h(x) 



(3.56) 



S(x) 



cr 



n 6 

^2 



which in turn defines the Bose invariant as 

b 2 2 

(3.57) j( x ) = _ + 

(j (j 



c 

The Schrodinger potential, given for z = J ^1 by 



1 - 



(3.58) 



t/(z) 



3 

4 rr 



4a 



1 



ft 6 

7^" 
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goes to oo as z — > oo and has a singularity at z = 0. The asymptotic behavior as z — ► implies 
that the Schrodinger operator has a self-adjoint closure with completely discrete spectrum. We 
have assumed |f (l — ^) < 1 which is always satisfied. Hence, the spectrum presents discrete 
eigenvalues of the form A„ = —bn for neN. The corresponding normalized cigenfunctions that 
satisfy the Schrodinger equation (2.341 are given in terms of the Laguerre polynomials 

a/a 2 



(3.59) 



<S>n{Z{x)) 



I>+f|) \2xJ 



-1/4 



2teY 



(H-i) 2b 



Theorem [16] yields the kernel for the semigroup generated by the operator C as a summation 
over Laguerre polynomials, which can be re-summed into (cf. (?), p. 705(7)): 

j>{T-t) \ Hft- 1 ) 



o(T-t)C 



(x,y) 



ye" 



exp 



— c I y + xe 



-b(T-t] 



' xye' 



-ft(T-t) 



(3.60) 

with c = c(T — t) = ||(1 — e _b( - T_ ^) _1 and in which one recognizes the transitional probability 
density of the CIR process. 

For the Laplace transform L, we specialize to the case where <p( x ) l& affile, i.e. 4>{x) = x. 
Therefore, (1.221 sets the shape of the polynomial Q(x,i9) to 

2d 



(3.61) 

The Bose invariant in this case is 

b 2 2d 



Q{x,ti)=S(x) + ^x 2 



(3.62) 



I{x) 



cr a*- 

whereas the Schrodinger potential, given by 



a b 

^2 



(3.63) 



U(z) = 



4o. 

7^2 



1 - 



1 , ab b 2 



2$a 2 
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corresponds to a discrete spectrum, following the same reasoning as for U (z). The solution 
is given in terms of the Laguerre polynomials 



(3.64) 



r(n+^) 



2x J 



1/4 



2bx\ 



a/a 2 



■2b 



with b — \Z2-da 2 + b 2 and for the corresponding eigenvalues A„ = —bn — -%•(£>—&). The Laplace 
transform L in ( 1.15 1 is easily integrated and yields a convergent series in terms of the Laguerre 
polynomials which re-sums to the famous formula (3.511 (cf. (?), p. 462(3) and p. 705(7)). □ 



20 



CLAUDIO ALBANESE AND STEPHAN LAWI 



4. Processes related to continuous orthogonal polynomials 



4.1. Proof of the Second Classification Theorem. We give a constructive proof of Theorem 
[5] independent of Theorem [2j In the following remark, we show that Theorem [5] can actually be 
regarded as a corollary of Theorem [2] 

The reducibility condition implies that the infinitesimal generator C must be of the form 

(4.1) C=^A {l) ^ + («-te)l 

for coefficients a € K and b > 0. The transitional probability density of the process, pT-t(x,y), 
satisfies the backward Kolmogorov equation: 



(4.2) 

with final time condition 
(4.3) 



dp , , , dp a 2 , , , d 2 p 



iim pr-t( x >y) = 5 ( x - y)- 



By the reducibility assumption, a general solution to this equation is given by the following 
eigenfunction expansion in terms of the orthogonal polynomials Q n (x; a, b) 



(4.4) 



P= ^ h n {t)Qn(x;a,b). 



n=0 



According to (1.271, the functions of time h n (t) satisfy the ordinary differential equations 
(4.5) h n + X n h n = 

which admits the general solution, for t <T, 



(4.6) 



h n (t) 



Z n 6 



A„(T-t) 



The coefficients z n are given by the final time condition (4.3 1: 

00 

(4-7) J2 > 



i(x;a,b) =S(x-y). 



Hence, multiplying on both sides by Q m (x; a,b) and the invariant measure p(dx), before inte- 
grating over the domain of x, leads to the result (1.29) by the orthogonality property (1.231. 

The Laplace transform Lx-tix, $) satisfies the backward Kolmogorov equation with differen- 
tial operator given by (1.131: 

(4.8) 



dL , , . dL a 2 . . . d 2 L n , /ST 
+ ( a _ bx ) _ + - A (x) = dcp(x)L 



dt 



dx 



with final time condition 
(4.9) 



lim L T -t(x, = q(x). 



For the sake of having a clearer constructive proof, consider the following ansatz for the Laplace 
transform 



(4.10) 

The previous equation then reads 



V(x)L. 



Y lh + (a ~ bx ^ V ' L + VL '1 + \\Mx)\{V"L + 2V'L' + VL") = HVL, 



(4.11) 

where the symbol ' denotes differentiation in the x- variable. The function V is chosen to satisfy 
(4.12) 



(a - bx)V + 2—A(x)V' = (a - bx)V 



LAPLACE TRANSFORMS FOR INTEGRALS OF MARKOV PROCESSES 



21 



for some parameters a £ M. and b > 0, which implies 

(a - by) - (a- by) 
a*\A(y)\ 



(4.13) F(a;)=exp 
The function 4>(x) is specified as follows: 
(4.14) 



dy 



V cr 2 V" 
H{x) = (a - bx) — + —A(x)—, 



which is equivalent to (1.301 with the assumption that C regroups all the constant terms. This 
choice yields the following partial differential equation for the function L: 



(4.15) 



dL ,_ 7 .dL 

m +i ~ a - bx) dx- 



° 2 A, ^ L n 



Following the same reasoning as for the transitional probability denstiy, a general solution to 
this equation is given by the following eigenfunction expansion in terms of the same orthogonal 
polynomials Q n (x;a,b), but with different coefficients, 



(4.16) 



L = 2_, h n (t)Q n (x; a, b). 



According to (1.271, the functions of time h n (t) satisfy the ordinary differential equations 
(4.17) h n + Xnh n = 

where A„ are the eigenvalues corresponding to Q n (x; a,b). The general solution, for t < T, is 



(4.18) 



h n (t) 



,A„(T-i) 



where the z n are constants. The latter equation for h n (t) wi th the explicit form of V(x) in (4.13 1 
and the expression of L in (4.161 gives the expected result (1.31 1 for the Laplace transform. 



(4.19) 



The coefficients z n are given by the final time condition (4.9 1 
z n Qn(%] a, b) = q(x) exp I - 



71=0 



(a - by) - (a - by) 



o*\A(y)\ 



dy 



Hence, multiplying on both sides by Q m (x) and the invariant measure p(dx), before integrating 
over the domain of x, leads to the final result (1.321 by the orthogonality property (1.231 and 
concludes the proof of Theorem [5j 



Remark 25. The reducibility condition implies that the operators C and C have the form 

.2 p,2 g 

+ (a — bx)- 



(4.20) 



C = 
£ = 



a " A, \ 9 



ff a, x 9 



dx 

2 - d 

+ (a — bx)- 



dx 



for possibly different coefficients a,a£M. and b, b > 0. Hence, setting 

a -by 



(4.21) 



R( x ) = and h(x) 



exp 



a 2 A(y) 



dy 



in the First Classification Theorem [1] yields the Second Classification The 
consider the gauge transformation T^ defined by 

a -by 



, [3| Moreover, 



(4.22) h(x) 
Then we have the following relation: 
(4.23) 



exp 



° 2 A{y) 



dy 



ThTr L — — A(x) ^-^r + (a — bx) - — h — „ . , > • 
h 2 y 'dx 2 y 'dx 2 &A(x) 



Notice further that the function V(x) of the previous proof is related to the latter two gauge 
transformations by V(x) — ■ 
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4.2. The Ornstein-Uhlenbeck process. 

Corollary 26. Assume that L and L are reducible in the sense of Definition [^] to Hermite 
polynomials. Assume also that the function 4>(x) = x. Then the transitional probability density 
of the Ornstein-Uhlenbeck process is given by: 



y, oo —bnfT— t) 

_ e -Mv^) 2 £ e —— Hn (z( X ))H n (z(y)) 

n=0 



and the Laplace transform is given by the following convergent series: 

oo 

(4.25) L T _ t (x,#) = e§*J2 e ~ >m{T ~ t)z "Mz( x ))> 

n=0 

where z(x) — \/ \(x — ?) and z(x) = \ \(x — f). The coefficients z n are given by: 



(4.26) ^ J q (x) e-i*H n (z(z))e-M*-%) dx. 



Proof. This case fits the classification scheme in Theorem [5] if one selects A(x) — 1. This choice 
implies 

(4.27) 4>{x) = C+^{(a- bxf -(a- bxf) . 

In order to reduce <f>(x) to an affine function, we choose b — b, a = a — ^f- and C = a 2 $°2 ■ The 
operator C, which has the form 

< 4 ' 28 > £ -TW + <- ta >S- 

has the Hermite polynomials H n (z{x)) as eigenfunctions with eigenvalues A„ = — bn. The 
invariant measure density and the normalization factor are respectively 



(4.29) w(x) = 



l_ e -^(*-i) 2 and dl = nir 



<7 2 ir 



which leads, by Theorem |5) to the formulation of the kernel of the semigroup generated by L as 
the convergent series in (4.241. The latter scries re-sums using Mehler's formula to give (3.281. 
By Theorem [5j the Laplace transform is given by (4.251. The coefficients z n are given by (4.261 
and lead to the result (3.451, once integrated and re-summed. □ 

4.3. The CIR process. 

Corollary 27. Assume that L and L are reducible in the sense of Definition to Laguerre 
polynomials. Assume also that the function 4>(x) = x. Then the transitional probability density 
of the CIR process is given by: 

,. nn * / S {2b\^ 2*^ _2 S ,^ n \ e -bn(T-t) 26 v r (*f-l),2& v 

(4.30) p T . t (x,y)=(^j y«> l e r(n+gg) Ln V^" 
and the Laplace transform is given by the following convergent series: 

(4.31) L T _ t (x,#) = e ^Ye-^ T ^hJ^ 1] (f 2 x). 

n=0 

The coefficients z n are given by: 
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Proof. This case also fits the classification scheme in Theorem [5] if one selects A(x) = x. This 
choice implies 



(4.33) 



c) = C + J- ((a - bx) - (a - te)) + ^ ((a - te) 2 - (a - to) 2 ) . 



2tfer 2 :r 



In order to reduce <j){x) to an affine function, we choose a = a, b = \j2 r da 2 + b 2 and C 



^ — 1). The operator £, which has the form 



(4.34) 



L = 



-Ydx^ + {a ~ hx) 



d_ 

dx' 



-(£-1), 26 



has the Laguerre polynomials (~ ^ a; ) 38 eigenfunctions with eigenvalues A„ = — bn if 

a > 0. The invariant measure density and the normalization factor are respectively 



(4.35) 



w(x) 



2b 



1 £2-t 
r 1 „ r J ' 



and tL = 



which leads, by Theorem [5| to the formulation of the kernel of the semigroup generated by C 



orem [5] 
ries inT| 



as the convergent series in ( |4.30 ) which re-sums to give (3.49). The Laplace transform is given 
by (4.311 and the coefficients z n are given by (4.321. The latter results lead to (3.511, once 
integrated and re-summed. □ 

4.4. The Jacobi process. 

Definition 28. The Jacobi polynomials P„ ' (x) are defined by the following Gaussian hyper- 
geometric function for 16 [—1,1]: 



(4.36) 



— n, n + a + (3 + 1 
a + 1 



1-x 



!, ? i-0,l,2, 



Definition 29. The Jacobi process is solution to the following equation: 
(4.37) dX t = (a- bX t )dt + ay/Xt(l-Xt)dWt 

with initial condition X t= o = xq G (0, 1). 

Corollary 30. Assume that L and C are reducible in the sense of Definition^to Jacobi poly- 
nomials. Then: 

p 2 -f3^ 



(4.38) 



4>(x) 



a 2 fa 2 -a 2 



x 1 — x 

2a 1 ^ 1 a 2 / u _\ i ^ i „ m j =r, 2a 



fora = %-l> -1,13 = £(6-a)-l > -1 and a = ft - 1 > -1,0 = £(&-a)-l > -1. 
TTie transitional probability density of the Jacobi process is given by: 

°° p -^n(n+a+[3+l)(T-t) 

(4.39) p T . t (x,y)=y a (l-yfJ2 ^ P^\l-2x)P^\l-2y) 



with normalization constant 
(4.40) dl 



d 2 



T(n + a + l)T(n + (3 + 1) 



(2n + a + (3 + l)T(n + a + (3 + l)n\ 

The Laplace transform is given by the following convergent series: 

00 

(4.41) L T _ t (x,#) = x^(l-x) i ^ i J2 e ~^ n{n+&+rJ+lKT ~ t)z nPi &j3) (l-2x)- 

71=0 

The coefficients z n are given by: 



(4.42) 



1 



d n Jo 



j{x) P^\l-2x) x^{l-x)^dx 
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with normalization constant 
(4.43) i 



r , _ r(n + a + l)r(n + /3 + l) 
~ (2n + a + + l)r(n + a + j3 + l)n! ' 



Proof. This case also fits the classification scheme in Theorem [5] if one selects A(x) — x(l — x). 
This choice implies 

(4.44) 4>(x) = C+^-( + + (a + (3) 2 - (a + f3) 2 



M \ x 1-x 

Since C is an arbitrary constant, we set it to C = (a + (3) 2 — (a + /3) 2 . The infinitesimal generator 
C, which has the form 

(4.45) C=^ X (l- X )^ 2+ (a- bx )l x , 

has the Jacobi polynomials Pi a '^\l — 2x) as eigenfunctions with eigenvalues A„ = —^-n(n + 
a + (3 + I) if a > ~ 1 and /3 > — 1. The invariant measure density and the normalization factor 
are respectively 

(4.46) «,(*)=*•(!-*)' and d\ = - r <" + , a ± ( " + P ±£ - , , 
v ; w v ' " (2n + a + /3+l)r(n + a + /5 + l)n! 

which concludes the proof by Theorem [5] □ 

4.5. The Dual Jacobi process. We introduce the dual Jacobi polynomials by applying the 
transformation x t— > Z{x) — x(2—x) to the Jacobi polynomials defined in the previous subsection. 

Definition 31. The dual Jacobi polynomials Dn (x) = P^ Q,/3 ' ) (1 — 2Z(x)) are defined as 
follows: 



(4.47) (x) = ^±^ 2 Ji " + a t ^ + 1 

n! \ a + 1 



x(2-x)\, n = 0, 1,2, 



They also satisfy an orthogonality relation (1.23) on (0, 1) with normalization constants and 
continuous measure density: 

T(n + a + l)T(n + /?+l) 



(4.48) w(x) = 2(x(2-x)) a (l-x) 2fi+1 , d 2 n 



(2n + a + j3+ l)T(n + a + (3 + l)n! ' 



They are solutions to the eigenvalue problem (1.271 with generator: 

a 2 , , d 2 2a - (2b - C)x(2 - x) d 
(4.49) C =-x(2-x) 1 ^ + 



dx 2 4(1 — x) dx 



and eigenvalues A„ = — ^~n(n + a + (3 + I), still conditioned to a = =§ — 1 > —1 and /3 = 
-% (6 — o) — 1 > — 1. Hence, we have the definition of the dual Jacobi process as follows: 

Definition 32. The dual Jacobi process is solution to the following equation: 

2a - (2b - C)X t (2 - X t ) a , — -, r 

(4-50) dX t = 1 4(1 Xt ) dt + ^ Xt{2 - Xt)dWt 

with initial condition X t= o — xq G (0, 1). 

We obtain yet another corollary to Theorem [5] 

Corollary 33. Assume that C and £ are reducible in the sense of Definition^to the dual Jacobi 
polynomials. Then: 

, , „ . a 2 fa 2 -a 2 '0 1 - 1 

(4.51) 4>{x) 



8^ \x(2-x) (1-x) 2 J ' 
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for a = |§ — 1 > — 1,/3 = -^(b — a) — 1 > —1. The transitional probability density of the dual 
Jacobi process is given by: 

(4.52) PT _ t (x,y) = 2(y(2-y)) a (1 - y) 2 ^ 1 £ ^—^(x)^'^). 

n=0 ™ 

TTie Laplace transform can be expressed by the following convergent series: 

oo 

(4.53) Lr-t(.x,0) = 2(x(2-x))^(l-xf~P ^ e x - {T -^z n D^^ {x) 
where the coefficients z n are given by: 

(4.54) z n = i / (a) (z(2 - z)) ^ (1 - a^+^da:. 

"n JO 

wrai/i normalization constant 

(4.55) d 2 - r(n + a + l)r(n + / 3+l) 



(2n + a + /3 + l)r(n + a + /? + l)n! 



Proof. The proof follows from Corollary 30 and the transformation x i— ► Z(x) — x(2 — x). □ 
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5. Processes related to discrete orthogonal polynomials 

The proof of Theorem [7] follows a very similar reasoning as the proof of Theorem [5j its 
continuous version. 

5.1. Proof of the Third Classification Theorem. The reducibility condition implies that 
the infinitesimal generator C must be of the form 

(5.1) C = -D{x)A 1 + (D{x) - B(x))V\. 

The transitional probability density pT-t(x, y) satisfies the backward Kolmogorov equation with 
generator given by (1.121: 

dt 



(5.2) ^-D(x)A 1 p+(D(x)-B(x))V 1 +P =0 



with final time condition 

(5.3) lira pT- t (x,v) = S(x-y). 

By the reducibility assumption, a general solution to this equation is given by the following 
cigcnfunction expansion in terms of the discrete orthogonal polynomials Q n (x) 

oo 

(5.4) P=X>»WOn(z)- 



According to (1.271, the functions of time h n (t) satisfy the ordinary differential equations 

(5.5) h n + X n h n = 
which admits the general solution, for t <T, 

(5.6) h n {t) =z n e A "( T -*). 



The coefficients z n are given by the final time condition (5.3 1: 

oo 

(5.7) ^ z n Q n (x) = S(x - y). 



n=0 

Hence, multiplying on both sides by Q m (%) and the weight w(x), before summing over the lattice 
Am, leads to the result ( 1.35[ ) by the orthogonality property (1.23). 

The Laplace transform LT-t(x, 1) satisfies this time a finite difference version of the Backward 
Kolmogorov equation 

clJ 

(5.8) — - D(x)A 1 L + (D(x) - B(x))V]_L = -d<\>L 

at x ' T 

with the same final time condition 

(5.9) lim £ T _t(x,i?) = q(x). 

Consider the ansatz for the Laplace transform 

(5.10) L = V{x)L. 
The latter finite difference equation reads 

dL V(x + 1) V(x — 1) 

(5.11) — - B(x) y > L(x + l)+ {B(x) + D(x)-^) L(x) - D{x) y > L(x - 1) = 0. 

V > ' B(x) + D(x) S ' 

B(x) D(x) 

B(x) and D(x) are defined such that they satisfy the relations 

V(x) B(x-1) D{x) 



(5.12) 



V{x-l) B(x-l) D{x) 
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which implies condition (1.371 and solves iteratively to give 

" D(k) 



(5.13) 



v{ X ) = n 



k=l 



D(kY 



The function 4>(x) is specified as follows: 



(5.14) 



1 



4>{x) = - (B(x) + D{x) - B(x) ~ D{x)) . 



(5.15) 



D(x)A 1 L + (D(x) - B(x))V\L = 0. 



But 4>(x), B(x), D(x) are all by definition independent of the parameter i9, so we are bound to 
set § = 1. This choice yields the following finite difference equation for the function L: 

dl 

In 

A general solution to this equation is given by the following eigenfunction expansion in terms of 
discrete orthogonal polynomials: 

N 

(5.16) L=Y^K{t)Q n {x). 

The functions of time h n (t) satisfy the ordinary differential equations 

(5.17) h n + X n h n = 
which admits the general solution 



(5.18) 



h n (t) 



,ln(T-t) 



where the z n are constants. The latter equation for h n (t ) with the explicit form of V{x) in (5.13 1 
and the expression of L in (5.161 yields the expression ( 1.38| fo r the Laplace transform. 
The coefficients z n are given by the final time condition (5.9 1: 

' : D{x) 



(5.19) 



n 



D(x) 

fc=l v ' n=0 



Finally multiplying on both sides by Q m (x) and the weight w(x), before summing over Ajy, gives 
the final result (1.391 by orthogonality of the polynomials and concludes the proof of Theorem 

13 

5.2. The Meixner process. The Meixner polynomials provide a discrete lattice approximation 
to the Laguerre polynomials. 

Definition 34. The Meixner polynomials are defined as follows in case x is integer: 

-n, —x . 1 





(5.20) 



M n (x;/3, c) = 2 Fi 



1 



!, n = 0,l,2,.. 



The Meixner polynomials satisfy an orthogonality relation with respect to the discrete measure 
supported on Z + . Namely, 



(5.21) 

where the weight is 
(5.22) 



^2 M m (x;l3,c) M n (x;(3,c) w(x) 



-"77.1 



x=0 



w(x) = — C . 

x\ 



The Meixner polynomials are solutions to the eigenvalue problem (1.271 with generator 

2 

(5.23) C = — x A 1 + (a- bx)V\, 
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for x £ Z + , with a, b > and eigenvalues A„ = — bn. The latter can be recast in the form ( 1.34 1 
using the functions 



(5.24) 



B{x) = - — c{x + (3) 
D(x) = — tt-X, 



for 



-/? and b = V(l — c). The parameters are conditioned to /3 > and < c < 1 which 



insures the Markov property. The Meixner process is the discrete Markov process generated by 
L. It is a discrete version of the CIR process. 

The following statement is a corollary to Theorem [7j 

Corollary 35. Assume that C and £ are reducible in the sense of Definition^ to the Meixner 
polynomials. Assume also that the function 4>{x) is given by 4>{x) = gx + C, 



(5.25) 



Q= y(c(e"-l) + (e-*-l)) , C = y/3c(e»-l) 



with the real parameter ip < — |lnc. Then the transitional probability density for the Meixner 
process is as follows: 



PT-t(x,y) = (i-c) p - — 



(5.26) 



■ 2 



F 1 \ 



(1 - ce (T-t)(c-l)-}x+y+f3 y \ 

-x,-y e (^)( C -D(i-c) 2 \ 
P c(l-e( T -*)( c - 1 )) 2 )' 



For q(x) = exp [u)(f){xy\, the Laplace transform is affine: 

(5.27) L T -t(x, 1,U) = e ™(T-t^)x+n(T-Uu) ^ 

The functions of time m(r; u>) and n(r; u>) are as follows: 



m(r; u>) = log 
n{r;uj) = -/31og | 
where c = ce 2tp . 



I — ce gu+ip _ ce(c-i)e-fT(l _ e gu-ipj J 

(c—l)e.~ v rH _ e QUJ-tp^ 



e 2 



cc(l-e*) 1 ~ CgW _ ,., 



1 - c 



Proof. The transitional probability density follows from equation ( 1.35 1 in the discrete classifica- 
tion theorem. The definition of 4>( x) suggests that we set B(x) = B(x)e v and D(x) = D(x)e~ v 
in order to satisfy condition (1.37) in Theorem u\ The generator defined by B(x) and D(x) has 
eigenfunctions M n (x;f3,c) with eigenvalues A„ = ^-e _v n(c— 1). Now from (1.391, we have 



(/3) n c"(l-c)' 9 



e^W"^ M n {x;P,c) 



( 1 


' "if 




VI- cee^-v J 




1 1 — ce^-v / 
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The Laplace transform, given by (1.38), is as follows: 



1 -c 

1 - ce eu} -f 

^ n! 

n=0 



/3 



ce 2 



1 gu — ip 

(c-l)e-V(T-t) 1 e 



1 - ce^-f 



M n (x;/3,c) 



/ i _ ce Qu+v _ (i _ eeu.-*.) e ^(s-i) e -"(T-t) \ * 

Note that the re-summation formula used to find the last two results is the generating function for 
the Meixner polynomials which can be found in (?). Also notice that M n (x; f3, c) = M x (n; (3, c) 
by definition. □ 

5.3. The Racah Process. 

Definition 36. The Racah polynomials R n (X(x)) :— R n (X(x); a,(3,~/,8) are defined as follows: 



(5.28) R n (\(x);a,P,'y,5) = 4 F 3 



-n,n + a + /3 + 1, —x, x + 7 + 8 + 1 



1 , »=0,1,2,...,JV 



where X(x) — x(x + 7 + 0" + 1) and either a = — TV — 1 or /3 + 8 = ~N — 1 or 7 = — iV — 1. 

The Racah polynomials satisfy an orthogonality relation with respect to the discrete measure 
supported on the set Ajv- Namely, 



(5.29) 

where the weight is 



n 



j;eA f 



(5.30) w(x) := w(x; a, {3, 7, 5) = 



and the normalization factor is 



(a + l) x {f3 + 8 + 1)^(7 + l) a (7 + S + 1) J; (( 7 + 8 + 3)72), 
(-a + 7 + 8 + l) a (-/3 + 7 + l)x((7 + <5 + 1)/2) B (<S + 



(5.31) 
with 



M 



(n + a + g + l) w (a + g - 7 + !)»(" - S + !)„(/? + l) n n! 
(a + /3 + 2) 2n (a + l)„(/3 + 5 + l)„(7 + l)„ 

r (-/3)jv(7 + <5 + 2) jv 



M = < 



(-(3 + j + l) N (S+l) N 
(-a + S) N (-f + 8 + 2) N 

(-a + -f + S + l) N {S+l) N 
(a + /3 + 2) N {-5) N 

(a-5+l) N (/3 + l) N 



if a = -N - 1 
if /3 + (5 = -AT - 1 
if 7 = -N - 1. 



The Racah polynomials are solutions to the eigenvalue problem (1.271 with generator (1.341 
given by the functions 

a 2 [x + a + \){x + (3 + 5 + l)(x + 7 + l)(x + 7 + 8 + 1) 



B{x) 
(5.32) D(at) 



2 (2a; + 7 + 8+ l)(2x + 7 + 5 + 2) 

cr 2 x(x - a + 7 + <5)(x — /3 + 7) (a: + 8) 
~2 (2x + 7 + (5)(2s + 7 + 5 + l) ' 



and eigenvalues X n = — %^n{n + a + (3 + 1). The Markov property is insured if B(x) < and 
D(x) < 0, Vx € Ajv- The process generated by the latter generator is called the Racah process. 
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Also define the corresponding functions 

cr 2 (x + a + l)(x+0 + 5 + l)(x + 7+ l)(x + 7 + 5 + 1) 



B(x) 

(5.33) D(x) 



2 (2x + j + 5+ l)(2x + 7 + 5 + 2) 

cr 2 x(x — a + 7 + S)(x — /3 + 7) (a; + 5) 
T (2x + 7 + ^)(2^ + 7 + ^+l) 
for a,/3,7,(5 £ R. 



Definition 37. TTie se£ of parameters {a, (3, 7, (5} toZZ &e called acceptable with respect to {a, /3, 7, 5} 
if it satisfies condition (1.37) in Theorem^ i.e. if B(x—l)D(x) = B(x—l)D(x), and ifB(x) < 
and D{x) < 0, Vz € A N . 

The following statement is a corollary to Theorem [7j 

Corollary 38. Assume that C and £ are reducible in the sense of Definition to Racah poly- 
nomials. Assume also that for an acceptable set of parameters {a, /3, 7, 5}: 

(5.34) 4>{x) = B(x) + D(x)-B(x)-D(x). 

Then the transitional probability density for the Racah process is given by 

N e -^n{n+a+l3+l)(T-t) 

(5.35) p T -t{x, y)=J2 J 2 R n (X(x))R n (X(y))w(y). 

while the Laplace transform is given by the following convergent series: 

(5.36) L T _ t (x,l) = UW) E e-^^^-^zMMx); a, 0, 7, 6). 

k=l ^ ' n=0 

The coefficients z n are as follows: 

(5.37) z n = i Y, ft <l{x)Rn{Xx); a, 0, 7, S)w(x) 

n xGAjv fe=l 

where d n = d n {a, 0, 7, 5), u)(a;) = w(x; a, 0, 7, 5) and A(x) = A(x; 7, 6). 



Proof. The restrictions imposed on the set of parameters {a, 0, 7, 5} ensures that condition ( 1.37 1 
in Theorem [7] is satisfied. This is a necessary condition. The rest of the corollary is a direct 
application of Theorem [7] □ 

5.4. The Dual Hahn Process. 

Definition 39. The dual Hahn polynomials R n (\(x)) := R n (X(x); 7, <5, iV) are defined as follows: 

(5.38) R n (\(xy, 7 ,6,N)= 3 F 2 (^ n >~*'*+^ 6+1 lj , n = 0, 1, 2, . . . , N 

where X(x) — x(x + 7 + S + 1) . 

For 7 > — 1 and S > — 1 or for 7 < — N and <5 < — N, the dual Hahn polynomials satisfy an 
orthogonality relation with respect to the discrete measure supported on the set A^r- Namely, 

(5.39) Rm(X(x))R n (X(x))w(x) = d 2 J nm 
where the weight is 

(2x + ^ + 6 + l)(j+l) x (-N) x N\ 
(-l)*(x + 7 + S + l)iv+i(<5 + l)xx\ 

and the normalization factor is 



(5.40) w(x) := w(x;j,6,N) 

and the normalization factor is 

(5.41) dl = 



1 



7 + n\ ( 8 + N — n 
N-n 
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The dual Hahn polynomials arc solutions to the eigenvalue problem (1.271 with generator (1.341 
given by the functions 

a 2 (x + 7 + l)(x + 7 + 5 + l)(N - x) 



(5.42) 



B(x) = 
D(x) = - 



2 (2x + -/ + 6+ l)(2x + j + 5 + 2) ' 
a 2 x(x + 7 + S + N + l)(x + S) 
Y (2x + j + 5)(2x + -f + 5+ 1)' 



and eigenvalues A n = — ^-n. The discrete Markov process generated by the latter generator is 
called the dual Hahn process. 

The following statement is another corollary to Theorem [7] 

Corollary 40. Assume that L and C are reducible in the sense of Definition to dual Hahn 
polynomials. Let 5 > 7 > — 1 or 5 < 7 < — N such that: 

^ r x(x + j + 8 + N + l) (x + ~/ + 5+ l)(N-x) 



(5.43) <j>(x) 



-(<5-7) 



(2.t + 7 + 5){2x + 7 + S + 1) O + 7 + 5 + 1)(2.t + 7 + <5 + 2) 



TTien £/ie transitional probability density for the dual Hahn process is given by 



(5.44) 



pr-t(x,y) = ^2 



e 2 



n(T-t) 



d 2 , 



-R n {\{x))R n {\{y))w{y). 



while the Laplace transform can be expressed as the following convergent series: 

N 



(5.45) 



L T _ t (x, 1) 



v ' > x n=0 



The coefficients z n are as follows: 



(5.46) 



(7 + 1), 



xEA N 



q(x)R n (X(x); 5, 7, N)w(x; S, 7, N) 



where d n = g?„((5, 7, TV). 

Proof. For 7 = 5 and (5 = 7, define the functions 

ct 2 (x + 7 + l)(x + 7 + (5 + 1)(AT - x) 
T (2x + 7 + ^+ l)(2x + 7 + ^ + 2) ' 
a 2 x(x + j + S + N + l)(x + 5) 



(5.47) 



D(x) = - 



2 (2x + 7 + (5)(2a; + 7 + <5+ 1) 
The corollary then follows from Theorem [7] 



□ 
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6. Limit relations 

Limit relations between orthogonal polynomials are well-known, see (?) and (?). In this 
section, we show that the transitional probability densities and Laplace transforms obtained 
from the various corollaries of Theorems [5] and [7] are similarly connected to each other. We start 
by rigorously stating the relation between the Jacobi process and its dual. 

Remark 41. The dual Jacobi process is obtained from the Jacobi process with the change of 
variable transformation 

(6.1) Z(x)=x(2-x) 

applied to the underlying process X t . The same holds for its transitional probability density and 
Laplace transform. 

For the discrete X t process, the limit relation between the Racah process and the dual Hahn 
process is not as obvious. 

Remark 42. The Racah process converges to the dual Hahn process in three different ways, 
corresponding to the three families of Racah polynomials: 

(1) a = —N— 1 with the acceptable set of parameters {—N— l,/3 + <5 — 7, £,7}, conditioned 
to either 

' (3>j + N ( p>-5-l 

<5 > 7 > -1 ° r \<5<7< -TV 

in the limit as (3 — > 00. 

(2) f3 = —5 — N — 1 with the acceptable set {a, —7 — -/V — 1,(5,7}, conditioned to either 

a>-f + 8 + N ( a > -1 

5 > 7 > -1 ° r \<S<7< -N 

in the limit as a — > 00. 

(3) 7 = — N— 1 with the acceptable set {— a + S — N — 1, /3, — N — 1, 6}, conditioned to either 

a > -1 ( a < -N 

/3>-l or I (i>-5-\ 

5>2a + N + 1 \5<2a + N+l 

with first the mapping St— >6 + a + N+ l and then the limit (3 — > 00. The dual Hahn 
parameters are in this third case (a, 5) . 

The result extends to the transitional probability densities and Laplace transforms. 

The next proposition states that the Meixner process converges to the CIR process in the 
affine case. 

x 

Remark 43. Under the transformations g 1— * g(l — c) and x t— > in the Meixner process, 

the limit c — > 1 yields the CLR process with parameter a = ft — l. The Laplace transform is affine 
in this case. 

Proof. The proof follows from the limit relation 

limM n ;a + l,c = " v ; . 

\l~c' J L { n a) (0) 

Both the transitional probability density and the Laplace transform in the Meixner case then 
converge to the affine CIR case, since 



□ 
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The last proposition demonstrates the connection between a discrete and a continuous un- 
derlying process X t by showing that the dual Jacobi process is actually a limiting case of the 
Racah process. 

Remark 44. Consider the acceptable set of parameters {a, (3, 7, £} where 



(6.2) 7 

Assume furthermore the inequalities: 



-(3 - S - N - 1 
—N — 1. 



(6.3) 

Then, applying the transformation x 
dual Jacobi process in the limit N — > 



Q 

7 



> 
> 
> 



-1 

P > 
S. 



i— > xiV io i/ie Racah process such that x € [0, 1], yields the 
oo for the special case of a = a and P > \f3\. This result 



applies to both the transitional probability density and the Laplace transform. 

Proof. The inequalities in the assumption insure that the process X t satisfies the Markov prop- 
erty as both B(x) and D(x) are negative for x £ Ajv- Since on top of a — a, 7 = 7, S = S, we 
have —f) + 7 = (3 + 5 and (3 + 5 = — (3 + 7, it is immediate that {a, (3, 7, 6} is acceptable. The 
function 4>{x) reduces to 
(6.4) 

x(x - a + 7 + 5)(x + 8) (x + a + l)(x + 7 + l)(x + 7 + 5 + 1) 
_(2x + j + 5)(2x + j + 5 + 1) ~ (2x + 7 + S + l)(2a; + 7 + (5 + 2) 



or equivalcntly 

= 

(6.5) 



-(/3-/3) 



x{x-2N -a- P- p- 2)(x - TV - 
(2x-2N-P-P- 2){2x -2N- 



-P-P-I) 
(3-13-1) 



(x + a + l)(x- N)(x -2N-0-P-1) 



(2x-2N - (3- (3- l)(2x - 2N - (3 - (3) _ 

4>{x) is bounded from below, from the inequality (3 > (3 > — 1 and is even monotonously increas- 
ing. Thus the assumptions of Corollary [38] are satisfied. So we have defined a proper Racah 
process with a Laplace transform that can be expressed as a convergent series over the Racah 
polynomials. With the transformation x 1— > Nx, the function cf>(Nx) converges, in the limit 
N — ► 00, to its continuous counterpart of Corollary 33 for $ = 1. 

Moreover, the finite difference infinitesimal generator C in ( 1.34[ ) converges to the dual Jacobi 
diffusion generator in ( 4.49[ ), under the transformation x Nx: 

(6.6) -D(Nx)A h + (D(Nx) - B{Nx))V+ 

a 2 2(a + 1)(1 - xf - (2(3 + l)x(2 - x) d a 2 x{2 - x) d 2 



2 4(1-1) dx 2 4 9a; 2 

in the limit /i = 1/N -> 0. 

Furthermore, the Racah polynomials, solutions to the finite difference equation generated by 



(1.34 1, converge to the dual Jacobi polynomials up to some factor 
n,n + a + (3 + 1 



4-F3 
(6.7) 

as N — > 00, with Z(x) 



iVx,iVx- 27V - 1 - 
a + l,-/3- JV,-JV 

n, n + a + /? ■ 
a + 1 



(a + 1), 



since 



x(2 



— * 2F1 
x). We also have as N 



1 



Z(x) 



00, 
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Finally, from all the above limit relations, the Laplace transform for the integral of the Racah 



process, given by (5.36) with the appropriate assumptions, provides an extension with underlying 



process X t on the lattice to the Laplace transform (4.53 1 for the integral of the dual Jacobi process 



in the particular case a = a, f3 > |/?| and (3 > —1. □ 

7. Conclusion 

We have given a complete classification scheme for diffusion processes for which Laplace 
transforms for integrals of stochastic processes and transitional probability densities can be 
expressed as integrals of hypergeometric functions against the spectral measure for certain self- 
adjoint operators. The known models such as the Ornstein-Uhlenbeck process, the CIR process 
and the geometric Brownian motion fit into this classification scheme. We have also presented 
extensions to these models in the quadratic Ornstein-Uhlenbeck process and the Jacobi process. 
An extension of the framework towards finite-state Markov processes related to hypergeometric 
polynomials in the discrete series of the Askey classification tree has been derived. Finally, we 
have explicitly computed some limit relations between discrete and continuous processes. 
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